Geometric Scaling in a Symmetric Saturation Modelfl 



Stephane Munier 



Universitd di Firenze 
Dipartimento di Fisica, via G. Sansone 1 
50009 Sesto Fiorentino, Italy 
munier@f i . inf n . it 



Abstract 



We illustrate geometric scaling for the photon-proton cross section with a 
very simple saturation model. We describe the proton structure function F2 at 
small- X in a wide kinematical range with an elementary functional form and a 
small number of free parameters. We speculate that the symmetry between low 
and high Q 2 recently discovered in the data could be related to a well-known 
symmetry of the two-gluon-exchange dipole-dipole cross section. 

1 Introduction 

Golec-Biernat and Wiisthoff (GBW) have shown that a saturation model is able to 
describe the HERA data for the proton structure function F 2 (x, Q 2 ) at small values of 
Bjorken variable x. The agreement of their model with the data is quite good for all 
available values of Q 2 , including the lowest ones. 

The GBW model is based on a dipole picture of photon-proton interaction 0. 
Its main ingredient is that at high energy, a particular frame can be chosen in which 
the interaction factorizes in two processes well separated in time. First the photon 
fluctuates in a qq pair of given transverse size r 7 , which then scatters off the proton. 
The factorized dipole-proton cross section a dp contains the QCD evolution. It only 
depends on the size of the dipole r 7 and on the Bjorken variable x. The following 
parametrization was chosen by GBW: 



where it was assumed that the transverse momentum scale Q s had a power-like depen- 
dence on x: 




(1) 




(2) 
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This choice was motivated by the following considerations. When rQ s /2<C 1, the model 
reduces to colour transparency. When one approaches the region rQ s /2~l, the expo- 
nential in Eq. (|1|) takes care of resumming many gluon exchanges, in a Glauber-inspired 
way. Intuitively, this is what happens when the proton starts to look dark. In the re- 
gion rQ s /2>l, such a microscopic interpretation is no more valid but the exponential 
has the virtue to force the cross section to tend to a constant at large rQ s (which 
also means at small values of x because of the ^-dependence of Q s ), thus respecting 
the unitarity constraint through a strict compliance with the Froissart bound. These 
formulae involve 3 free parameters only, o"o, Xq and A. Once those are determined by a 
fit to the F 2 data, one can take advantage of the universality of the dipole cross section 
to extend the model to more exclusive processes. Diffractive structure functions || are 
predicted in a quite satisfactory way. Elastic electro- and photoproduction of vector 
mesons is also accounted for very successfully @]. Some attempts were recently made 
to predict the 7* — 7* cross section |J. 

However, several points remain unsatisfactory in this model, both on the theoretical 
and on the phenomenological side. First, the form of the dipole cross section is quite 
ad hoc. Then, we note that for the description of small-Q 2 data, an arbitrary quark 
mass of order 140 MeV has to be chosen. On the phenomenological side, the original 
saturation model fails at describing the logl/x-slope of F 2 at large Q 2 . A revisited 
version has just been released ||, which takes care of some collinear resummations. 
The new model is much more successful, but it seems to incorporate less saturation 
effects since the saturation scale Q s is shifted towards a lower value 0. This might 
indicate that either saturation effects are there but can be mimicked by a DGLAP 
evolution, or DGLAP effects are still predominant. In the latter case the success of the 
original dipole model could be attributed to collinear resummations effectively taken 
into account in the dipole-proton cross section, though this is not at all apparent. 

Our goal in this paper is to go back to the main known features of saturation and 
to incorporate them in a simple parametrization for the structure function. Our main 
result is that an exponentiated elementary dipole-dipole cross section (see Eqs. (p!4] , |T5|) 



and Figs. 2,3,4) fits remarkably well the HERA data. We give some motivation and 
interpretation for this parametrization in Sec. |2| and we present the comparison to the 
data in Sec. |3|. 

2 Formulation 

Following GBW, we choose the so-called dipole-frame, in which all the QCD evolution 
is incorporated in the proton. The photon has just time to fluctuate in a qq pair 
which subsequently interacts with a highly evolved proton. We describe qualitatively 



the energy evolution of the proton in subsection (|2.1| ) and the scattering in Sees. (^2| 
and (|2.3|). The model is summarized on Fig.l. 



2 



2.1 The saturated proton 

When the center-of-mass energy increases (or equivalently when the Bjorken variable 
x decreases at fixed Q 2 ), more and more quantum fluctuations are revealed in the 
proton. Indeed, x coincides with the fraction of longitudinal momenta (with respect to 
the proton momentum), of the probed partons. The lifetime of a partonic fluctuation 
is proportional to this quantity. Thus going to smaller x means probing shorter time 
intervals, and thus becoming sensitive to more fluctuations. 

The rise of the parton densities with 1/x is predicted in QCD for a heavy onium 
||: the density n(x, k) of gluons of longitudinal momentum fraction x and transverse 
momentum k is found to obey a linear integro-differential equation, the BFKL equation 



H 10 1 . Its solution shows that the gluon density increases like a power of x: h(x, k) ~ 
x~ x . Appart from a diffusion term which would appear in a more refined treatment, 
this rise is independent of k. This model is not correct for very small x because the 
gluons interact with each other and thus each energy level can only accomodate a finite 
number of them: the density of partons of given transverse momentum should saturate 
at some point. This feature has been taken into account by a modification of the BFKL 
equation. This modification consists in supplementing it with a nonlinear term: 

dn(x,k) a s N c T ^ _ f __ ;x a s N c _ 2 



K-n{x,k) —h(x,k), (3) 

alog(l/x) 7T 7T 

where K is the linear BFKL kernel. This is the GLR equation ||11|| . Several groups 



have derived saturation equations within QCD (for a review, see Ref. [12|, |13[). Among 



them, Kovchegov was able to write a simple equation [14|] which reduces to Eq. @ 



when the partons are probed by a dipole of small size. All available approaches to 
saturation predict that the transverse momenta of the partons are on average shifted 
to a scale Q s called the saturation scale such that 

Q 2 s {x) = AV log(1/a;) . (4) 



This was confirmed by a numerical study presented in Ref. ||1 5|| . Corrections to for- 
mula (|j) have been computed recently, see Ref. |16| . The curve Q 2 = Q 2 (x) in the 
(x, Q 2 ) plane is called the critical line. The saturation scale Q s (x) is believed to reach 
1 GeV around x ~ 10~ 4 : this is what comes out of the GBW phenomenological ap- 
proach but a direct experimental evidence was also provided in Ref. |17| . 



It is convenient to pair the gluons in colour dipoles of given size r (which is possible 
in the large N c limit |§). The dipoles are then characterized by a density n(x,r). 
1/Q s (x) is their average size. It is also the mean distance between the center of neigh- 
bouring dipoles. We shall assume that the evolved proton is a collection of independent 
dipoles at the time of the interaction whose sizes are distributed around 1/Q s (x). We 
will approximate the distribution n(x, r) by a delta function of the mean value of the 
dipole size (r) = 1/Q S . Its complete evaluation is presently out of reach. Our confidence 



in this approximation relies on the observation made in Ref.|l8[ where it was shown 
that the data for a 7 * p for any x < 10~ 2 and Q 2 is function of the combined variable 
Q/Qs{x) only, to a very good accuracy. 
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To fix the normalization of n, we assume that the total area of these dipoles over the 
area of the proton is a constant ratio, of order 1. This is a kind of boundary condition 
on the critial line, in the spirit of Refs. [pTS|] and and yields the constraint 

fd 2 r Q 2 s (x) 

J ^n(x,r) = *±r 1 . (5) 

This assumption is supported locally by the saturation equations of the type @. How- 
ever, it also means implicitly that throughout the evolution, the area of the proton 
remains of order 1/A 2 . This can only result from non-perturbative confinement effects. 
It was indeed shown [[HJ that the Kovchegov equation, i.e. local saturation alone would 
lead to violations of unitarity due to the fast expansion of the nucleon radius. Note 
that the same hypothesis of a proton of fixed transverse size is contained in GBW 
model. 

The dipole distribution in the proton eventually reads 

n(x,r) = ^^-2nr-5(r-l/Q s (x)) . (6) 



2.2 The scattering 



In our model where the proton is represented by a density n of dipoles, the cross section 
for dipole-proton scattering reads 

a^(x,Q 2 ) = [ d 2 r, f dz\i> Q (r 7 ,z)\ 2 [ % n(x, r) • a dd (r„ r) , (7) 



o 



where ipQ is the photon wave function on a qq dipole state. 

a dd {r ryi r) is the dipole-dipole cross section computable in perturbative QCD. We 
only have to consider its lowest order expression, since all the QCD evolution is put 
in the dipole density n: there is no more room for radiative corrections in <j dd . It 
reads [f22j 

^(r 7 ,r) = 27ra s 2 r^l + log^ , (8) 

where r< = min(r 7 , r) and r> = max(r 7 ,r). The different factors that appear in 
this expression have a straightforward interpretation. This cross-section is the imag- 
inary part of the dipole-dipole elastic amplitude, which requires the exchange of at 
least two gluons paired in a colour singlet. These gluons must be able to resolve the 
smallest dipole, hence their transverse momentum is larger than l/r<. This is colour 
transparency. It accounts for the factor r<. The correcting logarithm, which becomes 
sizeable when the dipoles have very different sizes, comes from the monopolar Coulomb 
field, which results in a logarithm in 2 dimensions. Indeed, when r>3>r<, the smallest 
dipole sees the constituents of the largest one as two well-separated colour charges. 

In order to have an equal treatment for the photon and the proton and to get the 
simplest model, we treat the photon as a dipole of size 1/Q instead of taking account 
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of the complete probability distribution $(r 7 ) = J dzr 2 \ipQ(rj, z)\ 2 . This picture is 
relevant for a longitudinal photon whose wave function is sharply peaked (around 
logr 7 ~ log(l/Q s (a;))). It is less relevant for transverse photons since their wave 
function develops a plateau at large Q 2 ||23|| : this is due to the well-known fact that 



large size dipole configurations (aligned jet configurations) are always present in the 
transversely polarized photon. We shall however ignore this difficulty Indeed we 
believe it is not essential to our discussion since anyhow we donot pretend to be able 
to describe the very large Q 2 region. In this spirit, we approximate the distribution $ 
by a Dirac distribution: 

$(r 7 ) = N • 2vrr 7 • <5(r 7 - l/Q) . (9) 

N is a normalization factor that gives the rate of such a fluctuation 7* — > qq. 

We take a s fixed and at the scale 1 GeV. The approximation made here is that we 
stick to a kinematical region close to Q s where \\og(Q/Q s )\ <C \og(Q s / A). Injecting 
Eqs.(D, (|) and (§) into Eq.(|7|), we obtain: 

a^(x, Q 2 ) = N^^- ■ 2na 2 s (l + log for Q > Q s (x) (10) 

a^(x,g 2 ) = iV^-2vra^l+log^^) for Q < Q s (x) . (11) 

These formulae deserve an immediate interpretation. When Q > Q s (x) (first for- 
mula), we are in the usual picture of DIS where a photon of small size probes the 
parton (dipole) content of the proton. The flux factor is N, the target density is 
Q 2 (x)/A 2 , and the elementary cross section is 27ia 2 (l+\og(Q/Q s (x))/Q 2 . This is the 
hard Pomeron regime: the cross section grows like the exponential of the rapidity, due 
to the multiplication of dipoles in the proton. 

When Q < Q s (x) instead, the dipoles of the proton are smaller than the dipole of 
the photon. The use of a perturbative elementary dipole-dipole cross section is justified 
if x is very small so that the scale of the coupling constant Q s (x) is large enough. How 
large has always been subject to debate. We will assume that Q s ~ 1 GeV is fine. 
The proton looks like a collection of "small" size probes. The picture is inverted with 
respect to the previous case: 1/Q s (x), which is both the mean dipole size and their 
separation, fixes the resolution. The flux factor is Q 2 (x)/A 2 , the target density is N 
and the elementary dipole-dipole cross section is 2ita 2 (l + \og{Q s (x)/Q))/Q 2 (x). The 
growth of the parton densities with the energy is compensated by the photon-parton 
cross section, which falls with the energy. We are in the soft Pomeron regime, in 
which the 7* — p cross section grows only like a small power of the energy, (W 2 ) om . 
We note that in this region, the precise relationship between the photon virtuality 
Q 2 and the dipole size is not very relevant: corrections to it come as an additional 
energy-independent log term. 

An important propery of this parametrization is first that it only depends on 
Q/Q s (x) and second that the quantity Q /Q s {x) a 1 * v {Q /Q s (x)) is exactly symmetric 
under the exchange Q Q s (x). These two features have been observed in the HERA 
data P|. 
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We would now like to refine the model by including multiple gluon exchanges, which 
as we will argue, are important when Q ~ Q s ( x )- Throughout we will preserve the 
symmetry just pointed out. 



2.3 Accounting for multiple interactions 

Usually multiple gluon exchanges are disfavoured by the smallness of the coupling 
constant and/or by powers of Q 2 s (x)/Q 2 . We can see this by computing the probability 
that a dipole going through the proton undergoes an interaction. If the dipoles are 
uniformly spread over the surface of the proton, this probability is 

p = ^( X ,Q 2 ).^^2.a 2 ^ 2 (12) 

as long as it is small enough compared to 1. Note that p coincides with the packing 
factor k usually defined in the context of saturation (see for example [p4j|). 

We first observe that parametrically 2na 2 ~ 1 for a s evaluated at the scale Q ~ Q s ~ 
1 GeV. Then when Q 2 (x) <Q 2 , p is much less than 1: multiple gluon exchanges are 
disfavoured by a power of this small factor. When Q 2 ~Q 2 {x) however, p is of order 1 
and it is not justified to neglect multiple gluon exchanges. 

In a classical framework, these multiple interactions could be accounted for by 
modifying the interaction probability in the following way. We note that e~ p = S 2 is 
the probability of no interaction, where S is the S-matrix element for the dipole-proton 
interaction for a given impact parameter (it is constant over the area of the proton 
within our approximations). The dipole-proton total cross section is then 2-7r/A 2 - (1— S). 
Multiplying it by the rate N of splittings 7* — > qq, we obtain for the total photon- 
proton cross section 



a 7 p 



x,Q 2 ) = iV.|(l-5) = ^(l-e-/ 2 ) . (13) 



The exponential is reminiscent of a Glauber resummation which would be well justified 
for example in the case of a nuclear target. But we must keep in mind that the 
resummation is theoretically not under control for a proton target and there is up to 
now no systematic way to do it in QCD. 

We choose to keep the exponential form but we add a free parameter v > in its 
argument, so that we recover Eq.flTTD in the two-gluon exchange limit and the strength 
of the 4-gluon interaction be fixed by the data. We will check phenomenologically that 
the result is not very sensitive to u > . Absorbing all dimensionless normalizations in 
the parameters N and u > , we obtain the following formula: 

^^^{^(-.^m ( 1+log _«_))} . (14) 

Note that this formula is very close to the original GBW model (see Eq.(|l])). The main 
difference is the log present in the exponential, which will play an important role for 
the description of the data, as we shall see in the following section. 
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Let us now consider the case Q < Q s (x). We wish to preserve the symmetry of 
Q /Qs(x)a' y * p (Q/Q s (x)) for Q <->• Q s (x) at large \log(Q/Q a (x))\. Hence we take the 
following ansatz for the cross section in this regime: 

ff 7P (Sj g ) = JV ___| 1 _ e ^^ < _^ 1+log _jj| . (15 ) 

We have kept the same parameters except v > which is replaced by v K . Indeed, there is 
no reason why the resummation should be the same for Q > Q s (x) and for Q < Q s (x). 

Let us now try to interpret formula ([151) . The probability that a specific dipole 
in the proton has one interaction is p = {a 1 * v {x J Q 2 ) /n) ■ Q 2 /(ixN). By interpreting 
the proton as a single effective dipole, with flux Q 2 (x)/A 2 and the photon as a set 
of dipoles with density N, we are led to Eq.flTS]) by the same argument of a Glauber 
resummation of multiple gluon pair exchanges. Of course, this is only an a posteriori 
interpretation which has to be taken with great care. Indeed, it is usually argued that 
for Q < Q s (x), the dipoles in the proton are strongly correlated. 



3 A comparison to the data 

In this section, we compare Eq.([L||) (for Q > Q s (x)) and Eq.(|l5|) (for Q < Q s (x)) to 
the recent ZEUS data for the F2 structure function both in the high j2~51 and low pB 



Q 2 regime, and for x < 10~ 2 . We assume the parametrization (^) for the saturation 
scale Q s (x). We end up with 5 free parameters to be determined by the fit: Xq and A 
parametrize Q s (x), N is the global normalization, and v > and v < is the strength of the 
4 gluon coupling with respect to the 2 gluon coupling, for Q > Q s (x) and Q < Q B (x) 
respectively. 

The fit is performed with the 177 experimental points which remain after the se- 
lection x < 1CT 2 and Q 2 < 150 GeV . The x 2 is 1.09 per degree of freedom. It rises 
to 1.15 when we take all available Q 2 (185 data points). The parameters for Q s (x) 
are A = 0.35 and xq = 1.88 ■ 10~ 3 . The exponent for the rise of the parton densities 
A is consistent with the one found within other approaches. The saturation scale is 
1 GeV already for relatively high x (~ 2 • 10~ 3 ), but we cannot make a point of this 
fact since many neglected effects could shift the saturation scale by a constant factor 
(for instance the full photon wave function would have such an effect). The exponents 
za> and v < are very loosely determined by the data: we find v > = 0.8 and ^< = 0.6, 
with an error of about 50% and 100% respectively. If we fix them, to 1 for instance, 
the quality of the fit does not drop much (x 2 = 1.4/d.o.f.). 

The result of the fit is presented on Figs. 2 and 3. We have added some of the 
available data points which were not included in the fit and which sit at large Q 2 and 
large x. We see a good agreement over the whole Q 2 , x range. The x-slope seems 
to be reproduced fairly well even at large Q 2 . However, we slightly overshoot the 
normalization in this region. 

We have also shown, for comparison, the cross section without multiple scatterings, 
given by Eqs. ([X0|) and (|Tl"|). We observe that the data is well reproduced for Q<Q s (i) 
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and for Q Q s (x), but the resummations are needed around the transition region 
Q r ^Q s {x). We note that without the resummation of multiple gluon exchanges through 
the exponentiation, we would obtain a very poor global fit. However, we have a very 
good agreement with the data when we fit separately formulae (|10D and fllTD to the 
high and low Q 2 data, but with very different parameters x and A in the two regions. 

We plot the quantity Q /Q s {x)a' 1 *' p {x 1 Q 2 ) in Fig. 4. The symmetry under the ex- 
change of Q and Q s (%) is apparent. Our parametrization reproduces it exactly, by con- 
struction. We see that the Glauber resummations help to describe the early turnover 
seen near Q 2 /Q 2 (x) ~3. 

A final remark concerns the comparison with the Golec-Biernat and Wusthoff 
model. We have already noted that in the high-Q 2 domain, the main difference between 
our formula and the GBW model is the log in the exponent. We have checked 
that this term is crucial to fit the data. In the GBW model of Ref. 0, it is absent from 
the exponent, but there is an overall log factor coming from the photon wave function. 



4 Summary and discussion 

We can summarize our model as follows. We have assumed that in the dipole frame 
the energy evolution of the proton leads to the multiplication of the partons and conse- 
quently to the appearance of a transverse momentum scale Q s (x). When Q > Q s (x), we 
have the usual DIS picture, where the photon probes a set of independent partons. The 
rate of the growth of the parton densities has been parametrized by Q 2 (x)/A 2 . This 
procedure led us to formula ([H]). Then, we have taken seriously its extrapolation to 
small values of Q 2 , and we have observed that the symmetry of the dipole-dipole cross 
section by exchange of the dipole sizes is exactly seen in the data for 7*p total cross 
section. When Q~Q s (x), we have argued that the probability of multiple interactions 
between the photon and the proton becomes sizeable. A Glauber-like resummation 
was introduced for Q>Q s (x) and extrapolated by symmetry for Q<Q s (x). Although 
very simple and based on very qualitative ideas, the model agrees quantitatively with 
the recent high precision ZEUS data. 

This model provides an intuitive picture of geometric scaling both in the large 
Q (hard Pomeron, perturbative) and small Q 2 (soft Pomeron, non-perturbative but 
weakly coupled) regimes based on a symmetry between low Q 2 and high Q 2 recently 
observed in the photon-proton cross section. 

On the theoretical side, we may observe that our model could be related to unitary 
models based on Regge theory (see Ref.|27[]). In these models, a separation scale 
of 1 GeV -1 between small and large distance is assumed. This could correspond to 
our "saturation" radius 1/Q S which is of the same order. In our approach, the fan 
diagrams involving triple Pomeron vertices are hidden in the parametrization of the 
proton evolution (they are explicitly present in Kovchegov's approach for instance). 
The formula ([14]) may be seen as an eikonalization of fan diagrams. On the other 
hand, in a recent paper |T^] the dipole-proton scattering amplitude was derived in 
the context of linear BFKL supplemented by absorptive boundary conditions (see also 
The final formula looks very much like Eq. (|TO|) . It differs only by an anomalous 
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dimension which we donot take into account in our qualitative approach. The log 
has also presumably a different origin. However, this gives us confidence that the 
formula we propose is anyway well rooted in QCD in the restricted perturbative domain 
Q>Q s (x). 
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Figure captions 



Figure 1. A picture of the photon-proton interaction in the dipole frame 

In this frame, the momenta of the photon and of the proton are collinear. The 
photon splits in a qq pair of size l/Q. The QCD evolution leads to a set of dipoles of 
sizes 1/Q s (x) at the proton level at the time of interaction. These two systems scatter 
through the exchange of colour singlet pairs of gluons. The cross section is depicted 
for a. the exchange of two gluons, b. the exchange of four gluons. 



Figure 2. The F 2 structure function for medium Q 2 



Our parametrization of F 2 is shown (solid line) together with the latest published 
The fit was done only for x < 1CT 2 . 



ZEUS data 25 



The dashed line is the simple 
dipole cross section without taking account of many gluon exchanges (see formulae fllCf ) 
and (|Tl~D). The errors shown for the data are the quadratic sum of the statistic and 
systematic errors. 



Figure 3. The F 2 structure function for low Q 2 

The same in the low Q 2 region. The data points are from the ZEUS collaboration 
2l| . We only show the Q 2 -bins for which a significant number of data points are 
available. 



Figure 4. Geometric scaling 

The scaled photon-proton total cross-section Q/Q s (x)a ltp (x, Q 2 ) is plotted against 
the combined variable Q 2 /Q 2 (x). The data are from the ZEUS collaboration |21| , 
and are selected according to x<10~ 2 . 
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